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An approach is proposed to nuclear pairing at finite temperature and angular momen- 
, turn, which includes the effects of the quasiparticle-number fluctuation and dynamic 

coupling to pair vibrations within the self-consistent quasiparticle random-phase approx- 
imation. The numerical calculations of pairing gaps, total energies, and heat capacities 
are carried out within a doubly folded multilevel model as well as several realistic nuclei. 
^— ^ The results obtained show that, in the region of moderate and strong couplings, the 

sharp transition between the superconducting and normal phases is smoothed out, caus- 
ing a thermal pairing gap, which does not collapse at a critical temperature predicted 
by the conventional Bardeen-Cooper-Schrieffer's (BCS) theory, but has a tail extended 
to high temperatures. The theory also predicts the appearance of a thermally assisted 
pairing in hot rotating nuclei. 



1. Introduction 

The effect of temperature and angular momentum on pairing properties is an in- 
teresting subject in the study of nuclear structure. Because of its simplicity, the 
BCS theory is often used, which offers a good description of pairing correlation in 
the macroscopic systems such as metallic superconductors. It predicts a collapse of 
the pairing gap at T c , which signals the sharp supcrfiuid- normal (SN) phase transi- 
tion at finite temperature. The BCS theory, however, ignores quantal and thermal 
fluctuations, which are significant in finite small systems. Therefore, it needs to be 
corrected for the application to finite nuclei. Various theoretical approaches have 
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been proposed to study the effects of fluctuations on nuclear pairing l ^ l **! 4 ! . Their 
results show that, at zero angular momentum, thermal fluctuations smear out the 
sharp SN phase transition, resulting in a pairing gap, which does not collapse at 
finite temperature. In rotating nuclei, a phenomenon of temperature induced pair 
correlations, which reflects the strong fluctuations of the order parameter in small 
systems, has also been predicted The recent microscopic approach, called the 
modified BCS (MBCS) theory^ has shown, for the fist time, that the microscopic 
source causing the non-collapsing pairing gap is the quasiparticle-number fluctua- 
tion (QNF). 

Recently, we proposed the self-consistent quasiparticle random-phase approxi- 
mation (SCQRPA) 12, which includes the QNF as well as the quantal fluctuations 
due to dynamic coupling to pair vibrations. The purpose of present work is to 
extend this approach to finite temperature and finite angular momentum. 



2. Formalism 

The pairing Hamiltonian is considered, which describes a system of N particles 
interacting via a pairing force with the parameter G and rotating with angular 
velocity 7 and a fixed angular momentum projection M on the laboratory (or 
body) fixed z— axis: 

H = J2 e k (N k +N^ k )-Gj2Pl p k'-^N-lM , N ±k = a\ k a ±k , P\ = a\al k , 

(1) 

where a± k (a± k ) is the operator that creates (annihilates) a particle with angular 
momentum k, spin projection m k or —m k , and energy e k . For simplicity, the sub- 
scripts k label the single-particle states \k,m k ) with m k > 0, whereas — k denote 
the time-reversal states \k,— m k ). The particle number operator N is defined as 
ft = J2 k ( a l a k+ a - k a - k ), whereas M = J2 k m k (a\a k - a)_ k a^ k ) is the z-projection 
of total angular momentum. The variational procedure is applied to minimize the 
expectation value of this Hamiltonian in the grand canonical ensemble. The result 
yields the final equations for the pairing gap, particle number and total angular 
momentum, which include the effect of QNF in the form 

A k = A + SA k = Gj2^v k ,(V k ,) + G^\u k v k , (2) 
v { k) 

N = 2J2 Vl>P k ) + \{l-{V k ))\ , M = Y J ™ k {ni-n^) , (3) 

k L J k 

where the quasiparticle energy and renormalized single-particle energy el are 
given as 



E k = J(e' k -Gvl-Af+Al, (4) 
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4 = ^ + ^yE(^'-^')((4^) + (44')^) > (5) 

with (T> k ) = 1 — n k — n k , and A k = a\.a ! _ k . The expectation values (A\.Ak') and 
{A^^^) are evaluated by solving a set of coupled equations, which contain the 
SCQRPA X and Y amplitudes. The QNF is given as <SA" fe 2 = n+(l - n+) + n k (1 - 
n,7), where the quasiparticle occupation numbers n±fe are found from the integral 
equations 

7 ±H( e ^ + i)- 1 



r dw , (6) 



, [u-E k ± imk - M k (u)] 2 + fr^MP 

with the mass operators M k (ui) obtained by solving the set of equations for 
double-time quasiparticle Green's functions and those of a quasiparticle coupled 
with SCQRPA pair vibrations. The quasiparticle dampings are given as 7^(w) = 
SmfM^wifc)]. 

The proposed approach is called the FTBCS1+SCQRPA theory. Neglecting the 
coupling to SCQRPA, i.e. the factors (A f k A k ') and {A\A\,), it becomes the FTBCS1 
theory, which is different from the conventional FTBCS theory by the presence of 
the QNF. The violation of particle number at zero angular momentum is approxi- 
mately removed by applying the Lipkin-Nogami (LN) method. The corresponding 
approaches are called the FTLN1+SCQRPA and FTLN1. 



3. Results 

The numerical calculations are carried out within the il doubly degenerate equidis- 
tant model with the number £1 of levels equal to that of particles, N, as well as 
for 20 O, 44 Ca, 56 Fe, and 120 Sn. The results obtained show that, at zero angular 
momentum, under the effect of QNF within the FTBCS1 (FTLN1), the sharp SN 
phase transition predicted by the FTBCS theory is smoothed out. As the result, 
the pairing gap does not collapse at T = T c , but has a tail, which extends to high 
T . The dynamic coupling to the SCQRPA vibrations significantly improves the 
agreement with the exact results for the total energies and heat capacities obtained 
for N = 10 as well as those obtained 56 Fe within the finite-temperature quantum 
Monte Carlo method ^ [Figs. Q] (a) - Q] (c)]. However, for heavy nuclei such as 
120 Sn, the SCQRPA corrections are found to be negligible in comparison with the 
FTBCS1 (FTLN1) results. 

For 20 O and 44 Ca, the FTBCS1 pairing gaps, obtained at different M, decreases 
as T increases and do not collapses at high T. At M higher than the critical value 
M c , where the FTBCS gap for T = disappears, there appear thermally assisted 
pairing correlations, in which the FTBCS1 gap reappears at a given T\ > 0, and 
remains finite at T > T\ [Fig.[T](d)]. This phenomenon is caused by the QNF within 
the FTBCS1 theory At T = 0, the QNF is zero, so the FTBCS and FTBCS1 gaps 
are the same as functions of M (or 7), and both collapse at M = M c . However, 
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Fig. 1. Left panels: Pairing gaps (a), total energies (b), and heat capacities (c) obtained within the 
FTBCS (dotted lines), FTBCS1 (thin solid lines), FTLN1 (thin dashed lines), FTBCS1+SCQRPA 
(thick solid lines) and FTLN1+SCQRPA (thick dashed lines) for neutrons in 56 Fe. Boxes and 
crosses with error bars connected by dash-dotted lines are results of Ref. 8 . Right panels: pairing 
gaps as functions of T at different M (d), and as functions of M (e) and 7 (f) at various T obtained 
within the FTBCS1 theory for neutrons in 20 O. 

with increasing T, the FTBCS1 gaps, which are obtained at different T, collapse 
at M > M c , and remain finite even at very high T, whereas those given by the 
conventional FTBCS theory vanish at M > M c and T > T c [Figs. ED (e) and[T](f)]. 
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